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Abstract 

We discuss the principles to be used in the construction of discrete time 
classical and quantum mechanics as applied to point particle systems. In the 
classical theory this includes the concept of virtual path and the construction 
of system functions from classical Lagrangians, Cadzow 's variational principle 
applied to the action sum, Maeda-Noether and Logan invariants of the motion, 
elliptic and hyperbolic harmonic oscillator behaviour, gauge invariant electrody- 
namics and charge conservation, and the Grassmannian oscillator. First quan- 
tised discrete time mechanics is discussed via the concept of system amplitude, 
which permits the construction of all quantities of interest such as commutators 
and scattering amplitudes. We discuss stroboscopic quantum mechanics, or the 
construction of discrete time quantum theory from continuous time quantum 
theory and show how this works in detail for the free Newtonian particle. We 
conclude with an application of the Schwinger action principle to the important 
case of the quantised discrete time inhomogeneous oscillator. 

1 Introduction 

THERE are various circumstances in mechanics where it is convenient or necessary 
to replace the continuous time (temporal evolution) parameter with a discrete pa- 
rameter. Computer simulation of waves is an example where the configuration of a 
system at time t + T is calculated from a knowledge of its configuration at times t 
and t — T. There have been various attempts to construct classical and quantum 
mechanical theories based on this notion, such as the work of Cadzow |T[], Logan ||, 
Maeda || and Lee 0]. The work of Yamamoto et al. |5|, || and Klimek |7j indicates 
that the subject continues to receive attention. 

This paper considers the question: by which principles if any should continuous 
time mechanical theories be discretised, that is, turned into discrete time analogues? 
By discretisation we do not mean the numerical approximation of continuous time 
mechanics such as the work of Bender et al. [§]. Neither do we discretise space 



or the dynamical degrees of freedom. Our attention is fixed solely on replacing a 
continuous dynamical evolution parameter with a discrete parameter. In this and the 
following paper, Paper II on discrete time classical field theory ||, our interest is 
the construction of exact, self consistent discrete time mechanics with well specified 
principles, equations of motion and predictions. This is motivated by the notion 
that at some unimaginably small scale, time is really discrete. This has echoes in 
modern theories such as string theory and quantum gravity, where the Planck time 
of 10 -43 seconds sets a scale at which conventional notions of space and time break 
down. 

It could be argued that relativity requires a symmetrical treatment between time 
and space but this leads to the situation of a space-time lattice approach which has 
lost all relativistic symmetries and rotational invariance. We argue that relativity 
does distinguish between timelike and spacelike, and by discretising only time our 
approach reduces the break with relativity to a minimum. Lorentz covariance is 
broken in our approach to field theories, discussed in Paper II, but the residual 
Euclidean invariances permit the construct of particle like states. 

It may be felt objectionable that there is no natural concept of velocity in dis- 
crete time mechanics. It could be argued that this lack destroys our intuitive feeling 
for dynamics based on the notion of (say) a particle system evolving from an initial 
position and an initial velocity. The right way to see the situation is in terms of real 
numbers. In continuous time mechanics, we normally consider a particle as having 
an instantaneous position and an instantaneous velocity (we exclude Brownian par- 
ticle dynamics from our definition of continuous time mechanics). This information 
requires two real numbers for every degree of freedom. In discrete time mechanics, 
there is no natural concept of simultaneity analogous to this. What we mean by 
a "particle" is something with a position at time t and a position at time t — T, 
which also requires two real numbers for every degree of freedom. A particle here 
is more properly associated with the link between two successive points in discrete 
time, rather than those times separately. So ultimately, the only major difference in 
principle between continuous time and discrete time mechanics is the lack of the limit 
process T — > 0. 

One problem with discrete time mechanics is a lack of guiding principles at key 
places, which our series of papers attempts to address. For example, consider the dis- 
cretisation of a system with Lagrangian L = |mx 2 — V(x). The approach taken by 
most authors would be to replace the temporal derivatives by differences, symmetrise 
the potential in some way, derive the analogue of the Euler-Lagrange equation, and fi- 
nally evolve the system according to the resulting difference equation. Quantities such 
as the energy E = ^mx 2 + V(x) which are conserved in continuous time mechanics 
would be monitored by calculating the value Ed of the discretised Hamiltonian. 

It is more than likely however that a naive discretisation of the Hamiltonian would 
result in an expression E' D which is not exactly conserved. This has been discovered 
by many authors. It is a particular merit of Lee's approach [|J that an invariant 
analogous to the energy drops out of the formalism, but only at the expense of a 
dynamically evolving discrete time interval. 

It is somewhat surprising therefore that a computer simulation based on the above 
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principles should be judged good or bad according to how constant E' D remains. In 
the absence of any proper principle for the construction of invariants of the motion 
it should come as no surprise to find that occasionally a quantity such as E' D will 
vary enormously and unpredictably during the course of a simulation. This happens 
because there are really three systems being confused; 1) the original continuous time 
theory, 2) a discrete time system evolving exactly according to some well defined 
discretised Euler-Lagrange equation, and 3) some unknown discrete time system for 
which the naively discretised energy E' D would be an exact constant of the motion, 
but only for evolution under its own discretised equation of motion, which could be 
very different to the equation of motion for 2). On top of this there may be numerical 
uncertainties induced by the computer algorithms used. Seen in this light, it would 
seem a wise policy to discretise according to definite principles which would establish 
conserved quantities rigorously. The construction of invariants of the motion therefore 
becomes one of the principal objectives of discrete time mechanics. 

An important first step in the process of constructing rigorous discrete time me- 
chanics was the introduction of a discrete time action principle. This was done by 
Cadzow [|J, giving a discrete time analogue of the Euler-Lagrange equation. We 
shall call such an equation a Cadzow 's equation for the system. The construction of 
constants of the motion was considered by Maeda in the case of continuous sym- 
metries, whilst the construction of constants of the motion analogous to the energy 
had been considered earlier by Logan 0. 

Various features found in continuous time mechanics have discrete time analogues, 
including Noether's theorem, whilst certain other features either do not or cannot have 
discrete time analogues. An particular problem arises for example with Hamiltonian 
evolution and equations of motion derived using Poisson brackets. Not only is there 
no possibility here of an infinitesimal translation in time (which thereby renders the 
notion of a Hamiltonian problematical) but there is no natural concept of velocity 
as a limit either. This makes the standard definition of conjugate momentum as the 
partial derivative of the Lagrangian with respect to a velocity just as problematical. 
This has not prevented a number of authors from constructing discrete time analogues 
of Poisson brackets, however, with various degrees of success and utility, usually with 
the observation that the generator of time translations is not conserved. 

A feature of our approach is that we have found a clear principle for the definition 
of conjugate momentum in discrete time mechanics. It turns out not to be the 
partial derivative of the "Lagrangian" with respect to a difference in general, but 
does reduce to it in various important cases. In addition, we have avoided trying 
to construct Hamiltonians and equations of motion derived via Poisson brackets. In 
our formalism the Hamiltonian is displaced by a Logan invariant, if such a quantity 
can be found. Fortunately such an object does exist for the important case of the 
harmonic oscillator, which has ramifications in discrete time field theory discussed in 
the next paper in this series. 

The overall plan for this and subsequent papers is as follows. In this paper, 
Paper I, we restrict our attention to classical and quantum point particle dynamics, 
reserving classical field theory to Paper II, quantum field theory to Paper III, and 
quantum electrodynamics to Paper IV. 
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Topics covered in Paper I are as follows. First we introduce the central concept 
of system function. This replaces the Lagrangian as the key to the dynamics. With 
the system function we can calculate equations of motion, constructs invariants of 
the motion, and quantise the system. We give a prescription for constructing the 
system function from a given Lagrangian. We may use this prescription to embed 
into the system function symmetries such as gauge invariance and hence construct 
electrodynamics. Then we discuss the construction of invariants based on the work 
of Maeda, Noether, and Logan, and apply it to the harmonic oscillator, which we 
discuss in detail. This key system lies at the heart of particle field theory, discussed 
in the following papers, and it displays some important properties, such as a natural 
cutoff for particle energy for example. We also discuss particle electrodynamics and 
the Grassmannian oscillator. 

A major part of our programme is to develop discrete time quantum mechanics 
and so we conclude our paper with a discussion of the principles for first quantisa- 
tion. This includes the concept of system amplitude, the construction of unequal-time 
commutators, and compatible operators. We then discuss the construction of discrete 
time quantum mechanics from standard quantum mechanics via a stroboscopic ap- 
proach and give an explicit example. Finally we apply the Schwinger action principle 
to the discrete time inhomogeneous harmonic oscillator to construct the Feynman 
propagator for the oscillator, in anticipation of its use in field theory. 

2 Action integrals and action sums 

In continuous time mechanics Lagrangian dynamics is conventionally formulated via 
an action principle based on the action integral 



where t{ and tf are the initial and final times respectively along some given path T. In 
our version of discrete time mechanics we postulate that the dynamical variables q(t) 
are observed or sampled at a finite number of times t n , n — 0, 1, N, where to = ti 
and 1 7v = £/, such that the intervals t n+ i — t n are all equal to some fundamental 
interval T. For convenience we will write q n = q(t n ). 

It is possible to develop a theory where the time intervals vary dynamically along 
the path. Such a mechanics was considered by Lee [|J . The extension of our methods 
to that particular situation is left for a further article. 

In our formulation of discrete time mechanics we replace the action integral ([I]) 
by an action sum of the form 



where F n = F(c[ n , q n +i, n) will be referred to as the system function. The system 
function has the same central role in discrete time mechanics as the Lagrangian has in 
continuous time mechanics. With it we may construct the equations of motion, define 




(1) 



A N [T] = ]T F n , 



(2) 



n=0 
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conjugate momenta, construct constants of the motion and attempt to quantise the 
system. In principle we could consider higher order system functions which depended 
on (say) q n , q n+ i,..., q n+r , r > 2, but the case r = 1 represents the simplest possibility 
which could give rise to non-trivial dynamics and will be considered exclusively from 
now on. Such system functions are the discrete time analogues of Lagrangians of the 
canonical form L = L (q, q, t) . 

Another reason for considering only a second order formulation (r = 1) is its 
direct relationship to Hamilton's principal function, discussed presently. Cadzow |]J 
applied a variational principle to an action sum such as (Q) and derived the equation 
of motion 

A- Ipn-1 + F n\ = Q <n <N, (3) 

where the symbol = denotes an equality holding over a true or dynamical trajectory. 

We shall refer to ([|) as a Cadzow 's equation of motion for the system. We now discuss 
the interpretation of this equation. 

Suppose we have a continuous time action integral of the form (|IJ). First, partition 
the time interval [to^iv] into N equal subintervals. Then the action integral may be 
written as a sum of sub- integrals, i.e., 

4/[r]= E / dtL(q(t),q(t),t). (4) 

n=0 tn 

Now suppose that we fixed the co-ordinates q n at the various times to,ti, ...,ijv and 
then chose the path connecting each pair of points (q ra , q n+ i) to be the true or dy- 
namical path, that is, a solution to the Euler-Lagrange equations of motion for those 
boundary conditions. If this partially extremised path is denoted by T c then we may 
write 

N-1 

Af[T c ] = E ^ (5) 

71=0 

where S n = S(q n+ i, t n+ i, q n , t n ) is known as Hamilton's principal function, being just 
the integral of the Lagrangian along the true path from q n at time t n to q n +i at time 

tn+l- 

We recall now that the canonical momenta , p£+i at the end points q n , q n+ i 
may be obtained from Hamilton's principal function via the rule 

where the superscript (+) denote that the momentum at the initial time ti carries 
information forwards, whereas the superscript (— ) denotes that the momentum at the 
final time tf is influenced by earlier dynamics with respect to the temporal interval 
concerned. At this stage the action sum (|5]) has not been extremised fully, as the 
intermediate points q n , < n < N have been held fixed. 

Now suppose we went further and extremised ([5]) fully by variation of the pre- 
viously fixed intermediate co-ordinates q n ,n = 1,2, ...,N — 1. Then we would find 
that 

{S"- 1 + S n \ = 0, 0<n<N. (7) 
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This equation may be understood as the condition that the canonical momentum 
along the true path from to q/ is continuous, that is, p^ + ) = p£~). We notice 
immediately that (0) has the same formal structure as Cadzow's equation (|3]) provided 
we make the identification F n <-> S n . 

Another interpretation of Cadzow's equation is that it endows the action sum 
with the additivity property of action integrals, which satisfy the relations 

tl *2 *2 

J dtL + J dtL = J dtL, t < t! < h. (8) 

to tl to 

This property holds for all trajectories in continuous time mechanics, not just for the 
true or classical trajectory. In the case of system functions we may write 

F(q n _i, q n ) + F(q„, q n+1 ) = /(q„_i, q„+i) (9) 

for some function / of q ra _i and q ra+ i, because Cadzow's equation (Q) is equivalent 
to the statement that F n + F n ~ x is independent of q n along dynamical trajectories. 
However, unlike action integrals, this property will not hold off the true or classical 
trajectory in general. 



3 System functions from Lagrangians 

Two important ideas emerge from the similarity between (||) and (^): 

i) Although the concept of velocity as a limit does not occur in discrete time 
mechanics, we will define a unique discrete time momentum p n conjugate to q„ by 
the rule 

Pn = -^~F n . (10) 

»q„ 

This should be compared with the approach of Yamamoto et al. J5], || and most other 
workers, where the momentum is defined as a derivative of a discretised Lagrangian 
with respect to a difference. In our terms Cadzow's equation reduces simply to the 
statement that we may also calculate this momentum via the rule 

Pn = ^-F n '\ (11) 

ii) We will construct a system function F n from the temporal integral from t n to 
t n+ i of a continuous time Lagrangian, the question being which path to take. We 
cannot in general consider using the true continuous time path, as this is meaning- 
less in the context of discrete time mechanics and normally not known to us. For 
the particularly important case of the harmonic oscillator, however, we can evaluate 
Hamilton's principal function precisely and this provides us with an important check 
on our formalism. The path chosen in the construction of the system function will 
be referred to as a virtual path. 

It is possible to choose from a number of possible virtual paths, such as those 
inspired by q— deformed mechanics J7|. This does not alter any of the principles 
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we employ and simply changes the details of the system function used and hence 
the sort of invariants of the motion we can find. In this paper we are interested in 
treating time homogeneously, and so we choose a temporal lattice with a constant 
fundamental time interval T. Our proposed solution for the virtual path in point 
particle mechanics is to take the geodesic or shortest geometric path from q n to 
q n +i, the metric being normally the Euclidean one in physical space (not in co- 
ordinate space). This prescription will normally provides us with a unique system 
function from a given Lagrangian. Moreover, it should be applicable to configuration 
spaces with curvature and is a co-ordinate frame independent concept. It allows us 
to construct a gauge-invariant discrete time prescription for electrodynamics, with a 
suitable modification. In Paper II of this series we shall show that we can apply this 
prescription successfully to field theories also. There may be important cases where 
the chosen virtual path is not a linear interpolation. This occurs for charged fields 
in the next paper in the series. In such cases, additional requirements such as gauge 
invariance will influence the choice of virtual path. 

To illustrate the procedure, consider a non-relativistic particle with position vector 
x and Lagrangian L(x, x, t). Then the virtual path x n taken between x n and x n+i is 
given by 

x n = Ax„ +1 + Ax n , (12) 

where < A < 1 and A = 1 — A. With this choice of virtual path the time derivative 
become a difference operator. Specifically, we find 

~ d _ x n +i x n 

v n = -_x n - - , 

where we define 

t n = Xt n+ i + Xt n = t n + AT 
Then we construct the system function via the rule 

F n = T f 1 d\L(5t n ,v n ,t n ). (15) 
Jo 

The use of this integration does not imply that continuous time is regarded as mean- 
ingful in the context of discrete time mechanics. We are interested only in the results, 
not in the means of obtaining these results. A useful analogy is with the use of clas- 
sical mechanics to set up quantum mechanical models. Once we have found our 
quantum theory, we need no longer to regard the classical model which generated it 
as any more than some approximation useful in some circumstances. Our prescription 
allows us to embed into our system function fundamental properties such as gauge 
invariance and other symmetries of importance to physics. 

If the Lagrangian is a real analytic function of its arguments then we may make 
a Taylor expansion about x n and integrate term by term. This will be valid for 
Lagrangians which are polynomial functions of x and x. In such cases the system 
function F n would be given by the formal expression 

pn = T E , [ "\. L(x n ,v n ,t n ), (16) 
7 



(13) 



(14) 



where D n is the operator v n -^- + and v n and x n are considered independent at 
this stage. 

Some examples will illustrate the process. For a particle in one dimension with 
Lagrangian 

1 oo 

L = -mx 2 -Y,C r x T (17) 

^ r =0 

where the C r are constants, the system function is given formally by 

pn _ m ( x n+l ~ %n) _ r-p ~ X n~ ) (18) 

2T r t^ (r + l)(x n+1 -x n y 1 J 

For instance, the anharmonic oscillator Lagrangian 

L = —mx 2 muj 2 x 2 mXx 4 (19) 

2 2 4 K J 

gives the system function 

pn = m(x n+1 - x n f _ Tmu? (x 3 n+1 - xj) _ Tm\ (x b n+1 - x 5 n ) 
2T 6 (x n+ i - x n ) 20 (x n+1 - x n ) ' 

This differs from the anharmonic oscillator system function discussed in [ID|, which 
illustrates the general problem with discrete time mechanics. There may be many 
possible discretisation of a given continuous time system, all of which lead back to 
the continuous time theory when we take appropriate limits. The principle specified 
above gives us a unique discretisation (subject to choice of virtual path). 

For the coulombic potential problem in three spatial dimensions, the Lagrangian 

L= 2"*' i+ R (21) 

with virtual path (|12|) gives the system function 
m(x n+ i - x n ) 2 



2T 

_|_ 7^ ^ I x n+l'( x n+l — x n) + | x n+l| | x n+l — x n| 1 

| x n+l x n| ( x n'( x n+l x ri) ~l~ | x n|| x n+l x n| J 

This system function leads to Cadzow's equations of motion which preserve the dis- 
crete time analogue of orbital angular momentum. This system function is markedly 
different in form to the original coulombic Lagrangian ([H]) but if we consider tra- 
jectories for which we may write x ra = r n , x n+ i = r n + Tv n +0 (T 2 ) for each n, 
then ^ n 

lim | — 1 = imv n -v n + r ^ T , (23) 



t^o { T 

which corresponds with ( pi]) . However, it should be kept in mind that there will be 
many discrete time trajectories for which this limit cannot be taken. For example, 
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there may be trajectories where the particle repeatedly flips between two fixed posi- 
tions only. This may happen with the discrete time harmonic oscillator, for example, 
and no limit such as the one discussed above exists for such a trajectory Discrete 
time mechanics is inherently richer in its set of possible trajectories than continuous 
time mechanics. 

In general, Cadzow's equations lead to an implicit equation for x n+ \ involving x n 
and x n -i, although for certain systems such as the harmonic oscillator discussed below 
we may solve Cadzow's equation to find x n+ i explicitly. The situation is analogous 
to what happens in computer simulations of partial differential equations where not 
all equations give x n+ i explicitly. In such cases we must use numerical techniques to 
solve for the x n+ \ in the classical theory. It is a special feature of our approach that 
our equations of motion involve only x n -i, x n , and x n+ i, which is not always the case 
with finite difference schemes used to approximate differential equations. 



4 Invariants of the motion 

It is possible to find a discrete time analogue of Noether's theorem in the case of 
continuous symmetries along the lines considered by Maeda ||. We shall refer to 
constants of the motion found by this theorem as Maeda-Noether invariants . Consider 
a system function F n = F(q n , q n +i) which is invariant to some point transformation 
ch — ► q n = qn + Sq n . Then we may write 

Qpn Qpn 

= 5F n = -—8q n + - 5q n+1 

dq n oq n+1 

t " (24) 

using Cadzow's equation of motion. From this we deduce that the quantity C n = 
^-■8q n will be conserved along dynamical trajectories, that is, 

This construction does not allow us to construct an analogue of the Hamiltonian 
in the case of conserved systems because in our formulation we are not allowed to 
make infinitesimal jumps in time. 

Logan 0] gave a method for constructing constants of the motion which are not 
necessarily related to symmetries of the system function. Consider a point transfor- 
mation 

q« -> q„ = q«. + e u «, (26) 

where e is infinitesimal and u n is a function of q n and q n +i- Then 

f Qpn Qpn \ Qpn Qpn+l 

5Fn = e \ 7^- u « + ^ — u «+i = e 7^T' u " - e 7^ U «+!' ( 2? ) 

I oq n cq n+ i I c cq„ cq„ + i 



on the true trajectories. Now suppose that the transformation (|26|) is such that 5F n 
can be written in the form 5F n = ev n+ i—ev n , where v n = v (q„). Then we immediately 
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deduce that the quantity 

Qpn 

C n = —-u n + v n (28) 

is conserved over the classical trajectories. Such a constant of the motion will be 
referred to Logan invariant. 



5 The discrete time harmonic oscillator 

5.1 A Logan invariant for the harmonic oscillator 

The discrete time harmonic oscillator in its generic form is given by the quadratic 
system function 

F n = \a(x 2 n + x 2 n+l )- px n x n+u P>0, (29) 
which gives the Cadzow's equation of motion 

a 

x n+x = 2r]x n - x n -i, V=13- ( 30 ) 
6 P 

A Logan invariant of the motion is found to be 

C n = \fi {x 2 n + xl +1 J - ax n x n+ i. (31) 

5.2 Limiting behaviour 

In this subsection we show how to solve the equation of motion ( |30"D and determine 
the behaviour of the oscillator as the discrete time tends to infinity. First we define 
the variables 

a n = X n A* ^n+lj (32) 

which will become the analogues of annihilation and creation operators in quantum 
theory. The constants /i ± are chosen to satisfy the condition 

a n = A 4±a n— 1 ( 33 ) 

under the equation of motion ( ^0|) , which implies 

«n=(^) n 4- (34) 



Condition (|33D gives 



/i ± = rj ± Jr] 2 - 1. (35) 



We note that fi + fi = 1. The Logan invariant ([$!]) is given by 

C" = i/?a+a", (36) 

which is a constant of the motion by inspection and is a form of great value in discrete 
time field theory. 
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The complete solution to the problem is now readily obtained and given by 



. 1 ■;■ 



Xq 

-, ^1. (37) 



- fi-) 

For the case when rf < 1 we write ra = cos (8) and then we find 

_ sin (n6) x x - sin ((n - 1) 9) x 
c sin {&) 

whereas for rf > 1 we write rj = cosh (x) (assuming rj > 1), and then 

x = sinh (rax) - sinh ((n-l)x)x 
c sinh (x) 

and similarly for 77 < — 1. The crucial result is that bounded, elliptic behaviour occurs 
when rj 2 < 1 whereas unbounded (hyperbolic) behaviour occurs when rf > 1. This 
result gives a natural cutoff to particle momentum in field theory, as shown in Paper 
II. 

The readily solved case when rf = 1 corresponds to the free particle and will be 
referred to as the parabolic case. When rf > 1 the system will be said to be hyperbolic 
and elliptic when rf < 1. 



For the case 77 < 1 it is useful to define fj, — r} + iy/1 — rf and 



(40) 



the advantage being that these are constants of the motion, viz 

a n = a n -i, a* = a* n _ v (41) 



c c 



These are useful in the construction of particle states in quantum field theory because 
they correspond to annihilation and creation operators in the Schrodinger picture. 

5.3 The Newtonian oscillator 

Using the methods outlined in §3 the continuous time Lagrangian for the Newtonian 
harmonic oscillator 

^ 9 1 9 9 / \ 

L = -mi, — -muo x (42) 

gives the system function 

m(x n+1 - x n f TmuP 2 2 
F = ^ — (x n+1 +x n+1 x n + x n ), (43) 

which is equivalent to ( p9|) . The equation of motion is given by 

{%n+l 2x n -|- X n _i) 2 {%n+l ~\~ 4x n -|- X n — i) , . 

-UJ , (44 ) 



T 2 c 6 
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which is equivalent to fl30|) with the identification T ur = 6(1 — rj)/{2 + 77), which 
means 

6-2TV , v 

" = ITTrv ' < 45) 

Using the results of the previous subsection, we deduce that elliptic behaviour occurs 
only when the time T satisfies the condition 

< Tuj < 2^3. (46) 

An equivalent result is found in particle field theory, giving a natural cutoff for particle 
momentum. 



5.4 Harmonic recurrence 



We may understand the relationship between Hamilton's principal function for the 
interval [0, T] and the system function by explicitly evaluating the former for the 
continuous time harmonic oscillator Lagrangian We find 



2 sin(o;T) 

and comparing this with (^3|) we find 



'ra+l 



+ x 2 n ) cos(uT) - 2x n x n+ i 



(47) 



S n ( T ) = F n + Qf T 



(48) 

We expect a similar relation to exist in the general case, but different potentials will 
modify the precise details. 

There is an apparent problem with (57) whenever the time interval T satisfies the 
condition ujT = r7r , r = 1, 2, ... because the denominator sin (uT) vanishes at such 
times. This problem is an artefact of our representation of S n , because the definition 
of the principal function as a line integral over a finite contour of a bounded integrand 
means that S n cannot diverge. The resolution of this apparent paradox is that at the 
recurrence times T = rn/u the endpoints x n and x n+ i are no longer independent but 
are related by 

x n+1 = (-l) r x n . (49) 

The physical interpretation of recurrence is simple. The harmonic oscillator has 
a fundamental period P = 2h/uj , independent of the initial conditions. 

An important construction for the harmonic oscillator are the variables A n , A* n 
defined by 



A, 



A* 



xe 



sin (9) 

_ ie -ind 

sin (0) 



e- l6 x r 



9 = uT. 



(50) 



These are constants of the motion, i.e. A n = A n+ i and are independent of T. This 
means that recurrence must occur so as to cancel the zero of the denominator in fl4"7D 
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at the recurrence times. To see what happens explicitly, we may invert the equations 
(BOh to find 



. 1 .,- 



e me Al 



e i(n+i)e A * + e -i(n+i)e A 



so that at the recurrence times T = rir/u) we have 



(51) 



X ■)- 



[An + A r 



(52) 



from which we deduce (fj9|) . 

In terms of A n and A* n the principal function can be written as 



muj sin (8) 



e i(2n+l)8 A *2 _j_ e -i(2n+l)6»^2 



(53) 



a form which shows clearly that the principal function is not singular. Moreover, we 
see that at the recurrence times 



UJ 



0. 



1,2,... 



(54) 



It is a significant feature of the discrete time harmonic oscillator that it does 
not involve recurrence phenomena in this particular way, as no apparent singularities 
occur in the system function (^3|). This emphasises that discrete time mechanics is 
not equivalent to continuous time mechanics. 



6 Electrodynamics: test particles 

We now consider the case of electrically charged particles interacting with electro- 
magnetic fields. A more complete discussion of discrete time Maxwell's equations 
is given in Paper II of this series. Here we discuss only the case of test particles 
which are affected by external electric and magnetic fields but do not affect them. 
We shall find the Cadzow equation of motion for such particles and show that in our 
prescription electric charge is conserved. 

Consider a non-relativistic charged test particle of mass m in external electro- 
magnetic potentials. The continuous time Lagrangian for such a system is given 
by 

Lem = • x + gx • A(x, t) - g0(x, t), (55) 

where q is the charge of the particle. This Lagrangian is not gauge invariant but the 
equations of motion are gauge invariant, because under the gauge transformation 

<p' = + d t x, 
A^A' = A-Vx, (56) 

the action integral transforms according to the rule 

A if ^A! if = A if -[ qX ]%, (57) 
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that is, the change in the action integral occurs only at the end points. If this property 
is preserved in any discretised version of electrodynamics then the equations of motion 
should be gauge invariant. Our prescription for calculating the system function from 
the Lagrangian does indeed preserve this property and therefore our discrete time 
equations of motion are gauge invariant. 

The first step is to construct discrete time electromagnetic potentials. These are 
discussed in full detail in Paper 27, but the basic properties are the following. The 
magnetic vector potential A differs to the scalar potential <j) in that the former is 
defined at temporal lattice sites whereas the latter is defined on the links between 
these sites. If A n (x) is the value of the vector potential at time n at position x, and 
4> n (x) is the scalar potential on the link at position x from time n to time n + 1, then 
under a discrete time gauge transformation we have 

Xn+l ( X ) - Xn ( X ) 



</4( X ) = 0n( X )+ r 

A;(x) = A n (x) - V Xn (x) , (58) 

where \ n ( x ) * s the value of the gauge transformation function at time n and position 
x. The electric and magnetic fields are defined by 

E n (x) = -V0 w (x)- A ^W- A "W ; 



B n (x) = VxA n (x). (59) 

These are discrete time gauge invariant. By inspection the electric field is associated 
with temporal links whereas the magnetic field is associated with temporal sites. 

In order to apply our discretisation prescription to (^) we specify the virtual 
paths between times t n and t n+ i to be given by 

x n = Ax n+ i + Ax n , 
A n (x n ) = AA n+1 (x n ) + AA n (x n ) , 



(60) 



and then the system function is given by 

^ 2 rl _ rl 



F n = mlXn +i ^ + g(Xn+1 _ Xn) . / dX An (Xn) _ Tq I dX n ( Xn ) . (61) 



JO 



2T 

Under a gauge transformation we find 

F n ' = F n + q Xn (x n ) - q Xn+1 (x n+1 ) (62) 

N-l 

and so the action sum A = X) F n changes according to the rule 

71=0 

A N ' = A N + qXo -qXN, (63) 

in agreement with (|57D . Therefore we expect Cadzow's equations of motion obtained 
from ( 161]) to be gauge invariant. 
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In general, the integrals over the external electromagnetic potentials in ( |6lD give 
complicated equations of motion and we will normally have only an implicit equation 
for x n+1 , which however will be gauge invariant. We find 

m(x n+ i - 2x„ + x„_^ 



J>2 



q f dX {AE„ (x„) + AE„_! (x n _ x )} 
+g (Xn+1 T ~ Xn) x jf 1 dX {XXB n+1 (x n ) + A 2 B n (x n )} 



(64) 

In the limit T — > we recover the usual Lorentz force law 

mx = gE+gx x B (65) 

for those trajectories where the limit exits. 

The charge density p n (x) and current density j n (x) are defined by the following 
functional derivatives with respect to the electromagnetic potentials of the action sum 
A N : 

Pn (X) = y^TT^, jn (x) = I-^- A*. (66) 



We find 



p n (x) = q [ dX S 3 (x n - x) , 
Jo 

jn(x) = g (Xn+1 ~ Xn) fdXXPixn-x) 
1 Jo 



+g (Xn Xn ' l} rfA A 5 3 (x n _ x - x) , (67) 



T jo 

which satisfy the discrete time analogue of the equation of continuity 

Mx)-^_,(x) +v . n(x) = 



7 The Grassmannian oscillator 

We may apply our methods to the Grassmannian oscillator system, which serves as 
a prototype model for the Dirac equation studied in Paper II. Our model consists of 
one complex anticommuting degree of freedom ip with equation of motion 

i'ip = uip, i'ip = —u>ip*. (68) 

The Lagrangian giving these equations is 

L = ^i-ip*-ip — -i'ip ip ~ uip*tp. (69) 
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The equations ((68|) imply the harmonic oscillator equations of motion 
Now consider discretisation using the linear virtual paths 

i> = A^ n+ i + \lp n , ^* = X^n+l + ( 71 ) 

Using (|69|) we find the system function is 



F n = l -i 



2 



-wT- -, (72) 



which leads to the equation of motion 



n+l 



- c<j- -. (73) 



2T c 6 

and similarly for the complex conjugate. It seems not possible to use (|7"3|) to obtain 
the discretisation ( f44"D of the harmonic oscillator unless we change the virtual paths 
or renormalise the frequency to. However, we can readily show that ( [73]) implies 
harmonic oscillator behaviour by the following method. 
First, rewrite (f73f) as 

(_ 3 , + 6 ) ij n+1 + (3^ + 9) i, n _ x = -4^„, (74) 
where 9 = ujT. If we define 

v = v = (75) 

V9 + 9 2 V9 + 9 2 

and shift the degrees of freedom according to the rule 

= v n K (76) 
then the variables <p n satisfy the discrete time oscillator equation 

<p n+1 + n _! = 2r](f) n , (77) 

where n = V 2e . We find 



fi = V + i\/l-V 2 = r 1 » ( 78 ) 

9 + 9 2 



from which we deduce the upper limit ujT < v^3 for elliptic behaviour in the system. 
This is exactly one half of the upper limit found for the bosonic discrete time oscillator. 
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8 First quantisation 



We now discuss the possibility of quantising our classical discrete time mechanics. If 
we denote the process of quantisation by the symbol Q and the process of discretisa- 
tion by the symbol V then the question arises, do these processes commute, i.e. does 
QV = VQ. In other words, does it matter if we discretise the quantum theory of 
some system with classical Lagrangian rather than quantise the discretised version of 
the same classical Lagrangian? 

There is a priori no reason to expect these processes to commute. For one thing, 
we have not yet decided what Q might mean. Also, there are some aspects of V such 
as the lack of a Hamiltonian which makes a dynamical quantum theory of discrete 
time point particle mechanics problematical. Fortunately, there are some concepts 
from the conventional Q programme which are useful and appear to survive V. We 
shall comment on some of these aspects now and then consider the harmonic oscil- 
lator quantisation process QV in some detail in the following sections, reserving a 
discussion of the VQ process for Paper 77. 

In the following we discuss a system consisting of a point particle in one dimension. 
Generalisation to more degrees of freedom is straightforward and we shall use the 
Dirac bra-ket notation for convenience. 



8.1 Basics 

In our approach to quantisation, we shall follow all the standard principles of orthodox 
quantum mechanics in the main. This means we face the same issues of rigour 
and interpretation as orthodox quantum theory. We shall not comment on those in 
general. We discuss below those aspects where discretisation of time requires some 
additional emphasis or comment. 

Proposition 1: Physical states of a quantum system correspond one-to-one to 
rays in a separable Hilbert space 7i. A physical state vector \(f>) will be in general 
normalised to unity, viz 

(010) = 1- (79) 

Proposition 2: For each integer time n (or more accurately, at each co-ordinate 
time nT), H is spanned by an improper basis B n = {\x, n) : x G 3?} ; the elements of 
which satisfy the relation 

(x,n\y,n) = 5 (x — y) . (80) 
The resolution of the identity operator in H is 

which holds for each n. 

Given a physical state \ip) m.TL we may write for each n 

\ij) = J dxi> n (x)\x,n), (82) 
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where ip n (x) is the wavefunction at time n, with the property that 

dx\i, n {x)\ 2 = l 1 (83) 



assuming normalisation to unity. 

Remark 1: The Heisenberg picture is being used in the above. The time depen- 
dence of the basis sets B n allows us to use the Schrodinger picture, discussed below. 

Remark 2: All the usual principles of quantum mechanics concerning the in- 
terpretation of the states in the Hilbert space apply here. For example, a) the su- 
perposition principle and its interpretation according to standard quantum mechanics 
holds and b) given two physical states |0), then the inner product gives the 

conditional transition amplitude for the system to be found in state \<p), given it is in 
state 

Definition 1: An operator A diagonal with respect to B n is one which can be 
written in the form 

A= / dx \x, n)A{x 1 d x ){x, n\, (84) 



where the component operator A(x,d x ) is some differential operator of finite order. 
The action of such an operator on a typical state \ip) is given by 



\Aijj) = A\ip) = I dx \x, n)A(x, d x )ip n (x) (85) 
with matrix elements given by 

(0| Aij)) = I dx(p* (x) A(x, d x )ip n (x) . (86) 



Remark 3: The wavefunctions of the theory are elements of C 2 (9ft) , the space 
of square integrable functions on 9ft, and the operators (including the observables of 
the theory) which act on them are usually built up of functions of x and d x . If at 
a given time n the component operator of some observable diagonal with respect to 
B n happens to be represented by say a multiple of d x this carries no implication that 
the observable is related to a velocity (which would normally be implied in conventional 
wave mechanics, where the momentum operator is represented by —ihd x ). There is 
no concept of velocity in the normal sense in discrete time mechanics. 

Keeping in mind the caveats discussed by [|TJ[] concerning hermitian and adjoint 
operators, we define the hermitian conjugate or adjoint operator A + to have the 
property that 

(<t>\A^) = (^|i + 0>* (87) 

for a dense set of physical states. If A is diagonal with respect to B n then assuming 
we may represent A + in the form 



A + = J dx\x 1 n)A(x,d x )(x,n\ (88) 
for some operator component A(x, d x ) we find 

dx 0* (x) A (x, d x ) ip n (x) — dx A* (x, d x ) 0* (x) ijj n (x) . (89) 
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Assuming that we are permitted to integrate by parts, which will be the case for 
normalisable wavefunctions falling off at spatial infinity, we can readily understand 
the relationship between A (x, d x ) and A(x, d x ). Given the former we can always work 
out the latter by integration by parts and vice-versa. 

If A (x, d x ) = A(x, d x ) then A = A + and the operator is self-adjoint. Physical 
observables of the theory which are diagonal with respect to B n will be assumed to 
have this property. 

8.2 Dynamics 

The dynamical content of the theory is expressed in terms of unitary timestep oper- 
ators U n , one for each n. 

Proposition 3: For each n there is a unitary operator U n such that 

\x,n+l) = Ul\ x,n). (90) 

From this we deduce the relations 

(x,n + l\ = (x,n\U n , 

\x,n) = U n \x,n + 1), (91) 
(x, n\ = (x,n + l\Ul. 

Remark 4: The operator U n provides an isometry between B n to B n+1 , that is 

(x, n + l\y, n + 1) = (x, n\y, n) = 5 (x — y) . (92) 

Using the resolution of the identity ( |5T] ) we may represent the timestep operators 
in the non-diagonal form 

U n = J dx\x, n){x, n + 1|, — J dx\x,n + l)(x,n\. (93) 

We are now in a position to define the fundamental functions of the quantum 
theory, the system amplitudes U n (x,y), defined by 

U n (x,y) = (x,n + l\y,n) = (x,n\U n \y,n), 

U*\x,y) = (y,n\x,n + 1) = (y,n\Ul\x,n), (94) 

from which we arrive at the non-diagonal expressions 

U n = J dxdy\x,n)U n (x,y) (y,n\, 

Ul = J dxdy\x,n)U*(y,x) (y,n\. (95) 

Remark 5: The system amplitudes will not be differential operators of finite order 
in general but well-behaved complex valued functions of two real variables. Neither 
will they be singular distributions. They are similar in function to integrated Feynman 
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transition kernels encountered in the path integral formulation of standard quantum 
mechanics. 

The condition that the timestep operators are unitary, viz 

Until = In (96) 

leads to the closure condition 

JdyU n {x,y)U*(z,y) = 6(x-z) (97) 

on the system amplitudes. 

Definition 2: A system for which the system amplitudes are independent of time, 
such that we may write 

U n (x,y) = U(x,y) (98) 

for some U (x, y) and for all n, will be said to be autonomous. 

Definition 3: An autonomous system for which the system amplitude U (x, y) 
carries the symmetry 

U{x,y) = U(y,x) (99) 

will be said to be time-reversal invariant. 

Remark 6: Most of the system amplitudes of interest to us will be autonomous 
and time-reversal invariant. This will be so whenever we construct system amplitudes 
from system functions which have been obtained from conventional time-translation 
invariant and time-reversal invariant Lagrangians using the virtual path approach 
discussed above. 



8.3 The Schrodinger picture 

The Heisenberg picture description of the physical states used so far means that we 
may write 

\ip) = J dxip n+1 (x) \x, n + 1) = J dxift n (x)\x,n), (100) 
from which we deduce 

^ n+ i(x) = J dyU n (x,y)ip n (y). (101) 

From this we see that the system amplitudes play a role analogous to finite time scat- 
tering kernels in conventional quantum mechanics. Equation ( |101| ) is about the closest 
we come in this theory to something analogous to a time dependent Schrodinger wave 
equation. 

We may set up a formal description in the Schrodinger picture as follows. Given 
a Heisenberg picture state \ip) and a knowledge of the component functions ip n (x), 
define the sequence of states 

|V>„,J = / dx?p n (x)\x,m), (102) 
for some chosen time m. Then if \ip mm ) = 1^) we find 

Um\^ m , m ) = IVWl,m)- (103) 

It is straightforward to extend this to jumps over more than one time interval. This 
establishes the Schrodinger picture in this theory. 
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8.4 Position eigenstates 



Given an improper basis B 
tion operator 

which has the property 

Remark 7: For convenience we shall follow the traditional abuse of notation and 
use the symbol x both for the position operator and for a particular eigenvalue of that 
operator. It will be clear normally from the context what is meant whenever a clash 
of notation occurs. 

The position operators have the merit of being diagonal with respect to the appro- 
priate basis, that 1 S • X IS diagonal with respect to B n . The position operators are not 
necessarily diagonal with respect to bases at other values of n. From the condition 

x n+ i = Ulx n U n , (106) 

we find 

x n +i = J dxdydz \x,n)U*{y,x)yU n {y,z)(z,n\, (107) 

which is self-adjoint but not necessarily diagonal with respect to B n . 

Remark 8: It is possible for x n+ i to reduce to diagonal form with respect to the 
B n basis but this depends on the details of the system amplitudes. 

From the above we arrive at the fundamental expression for the commutators: 

[x n+1 ,x n ] = / dxdydz \x,n)U*(y,x) y (z - x) U n (y, z) (z,n\. (108) 



e {\x, n) : x G 3?} we may construct a self-adjoint posi- 
x n = J dx \x, n)x(x, n|, (104) 

x n \x,n) = x\x,n). (105) 



8.5 Normal co-ordinate systems 

In this subsection we discuss the quantisation of a large family of systems for which 
the following property holds: 

Definition 4: Co-ordinates for a system for which the right hand side of the 
commutator ( |108| ) is a multiple of the identity operator for each value of n will be 
called normal. 

Remark 9: We shall see below that the co-ordinates for the important system 
equivalent to the harmonic oscillator are normal. 

A class of system amplitudes for which the co-ordinates are normal may be con- 
structed from autonomous, time-reversal invariant system functions of the form 

F(x n , x n+ i) = -f3x n x n+ i + \W{x n ) + -W {x n+ i) , (109) 
where (3 is a non-zero constant. The Cadzow's equation for this system is 

X n+1 = f3' l W {X n ) - (110) 

which has the merit of giving x n+ i explicitly in terms of x n and x n -\. 
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Now define the system amplitude to be given by 

U n (x,y) = ke iF ^\ (111) 
where k is some constant. Then from the unitarity condition fl9~T|) we find 



\k\ 2 = —- (112) 
1 1 2nh K ' 



From (|108|) we find 

so that the above co-ordinates for this system are normal. Moreover, with the mo- 
mentum p n conjugate to x n defined by (10) , we recover the conventional commutator 

[p n ,x n ] = -%h. (114) 

This result is not expected to hold for systems which are not normal. 
The operator equations of motion are found to be 

x n + 1 = P~ 1 W' n -x n - X , (115) 

where W' is the diagonal operator 

r 4 f dw (x) 1 ; , , 

W' n = I dx \x,n) I ^ ' \ (x,n\. (116) 



dx 

From this we obtain the discrete time version of Ehrenfest's theorem; i.e. 

(x n+1 )=p- 1 (W; i ) -<£„_!> (117) 
for expectation values over a physical state. 

8.6 Compatible operators 

In our theory the quantum dynamics is completely determined by the system ampli- 
tude. Suppose now that the system is autonomous and time-reversal invariant. This 
means that for each time n we may write 

U n (x,y) = U(x,y), (118) 

where U (x, y) is independent of n. For such a system there may be constants of 
the motion comparable to the Maeda-Noether and Logan invariants discussed in the 
classical theory. 

Consider an operator C diagonal with respect to B n , viz. 

C= / dx \x, n)C{x 1 d x )(x, n\, (119) 
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where C(x,d x ) is some differential operator of finite order. Matrix elements of the 
commutator of C with U are given by 

(0| [C,U] |0> = J dxdyf n (x){C(x,d x )U(x,y)-C*(y,d y )U(x,y)}ij n (y), (120) 

where \ip) and \<f>) are arbitrary physical states. From this we arrive at the result 

Theorem: A diagonal operator commutes with the timestep operator of an au- 
tonomous system if 

C(x, d x )U (x, y) = C*(y, d y )U (x, y) . (121) 

Definition 5: A diagonal operator which commutes with the timestep operator of 
an autonomous system will be said to be compatible (with the timestep operator). 

Remark 10: It is not necessary for a diagonal operator C to be self-adjoint for 
it to be compatible with the timestep operator. 

Remark 11: From the above we deduce that compatible operators are invariants 
of the motion. To be explicit, consider a state \ip) which is an eigenstate of the 
diagonal operator C with eigenvalue c, i.e. 

C\fl;) = c#>. (122) 

Then we can show 

C{x,d x )^ n {x) = c^ n {x) (123) 

and 

C(x, d x )iJ n+1 (x) = cip n+1 (x) . (124) 

Remark 12: Given the system amplitude U (x,y) it may be very hard or perhaps 
even impossible to find any compatible operators in closed form. It may be necessary 
to approximate such an operator via a perturbative expansion, for example. This is 
the quantum theory analogue of the problem of finding invariants of the motion for a 
classical discrete time theory given some system function. 

Remark 13: Discrete time and continuous time quantum mechanics pose dual 
problems in the following sense. In continuous time quantum mechanics we are nor- 
mally given a Hamiltonian and the problem is to construct the time evolution opera- 
tor. For a time independent Hamiltonian a complete solution would require us to find 
all the eigenvalues E a and eigenstates \E a ) of H and then use them in the formal 
solution Ut = exp {-iHt/h} to write 

U t = Y,\E a )e-* E ^ h (E a \. (125) 

a 

This is in general a formidable problem. In discrete time quantum mechanics the 
situation is the other way around. Given a system amplitude, the problem is to find 
the compatible operators, if any. An important system where answers can be found 
to all of these questions in both approaches is the discrete time harmonic oscillator 
discussed in 510. 
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9 Stroboscopic construction 



In principle it should always be possible to construct examples of discrete time quan- 
tum systems by integrating the equation 

ihd t U{t) = HU (t) (126) 

for the evolution operator U(t) in continuous time quantum mechanics, given the 
Hamiltonian H. The boundary condition 

lim U(t) = i H (127) 

ensures a unique solution. From this point of view the discussion outlined above rep- 
resents a stroboscopic approach, where the state vectors evolve continuously but are 
looked at periodically in a non-destructive (mathematical) sense. This examination 
of the state vector is not the same as an observation collapsing the wavefunction. 
For autonomous systems a formal solution to (|126|) is 



U(t) = exp (iHt/h) . (128) 

In this approach the transition amplitude U(x,y;t) = (x,t\y,Q) corresponds to our 
system amplitude when t = T and may be evaluated in a number of ways. For 
example the Feynman path integral method gives the formula 

(x,t\y,0) ~ J [dz]exp^J dt'L(z,z,t')\ , t > 0, (129) 

where L is the Lagrangian, such amplitudes being functions of t and the end-points 
x and y. The standard approach to the evaluation of such integrals is, rather inter- 
estingly, based on the discretisation of time. The time interval [0,t] is partitioned 
into a finite number iV of equal steps, the integrand in the exponential is approxi- 
mated suitably (by what amounts to choosing a virtual path in our approach), the 
N integrals evaluated, and then the limit N — > oo taken. 

The relationship between this approach and our discrete time formalism should 
now be clear, the basic difference being that we do not take the limit N — > oo. In a 
number of situations our system amplitude will actually take the form 

U (x, y) = (x, T\y, 0) ~ exp {iF (x, y) /%} , (130) 

and then for N timesteps the transition amplitude (x, NT\y, 0) becomes 

(x, NT\y, 0) ~ J ... J dx\dx2-..dxN-i exp jzA^/ftj , (131) 

which emphasises the relationship further. 

An important point which could be confusing is that the system function in (|130|) 
does not correspond to a Logan constant for those systems such as the harmonic 
oscillator where such an invariant can be found. System functions in general are not 
expected to be invariants of the motion. 
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An alternative method of constructing the transition amplitudes is to find the 
Green's functions for the system. If the transition amplitude U(x,y;t) satisfies the 
homogeneous Schrodinger equation 

(ihdt-l£)U(x,y,t) = (132) 

with the boundary condition 

hm[/(x,y,t) = <5 3 (x), (133) 

then the retarded and advanced Green's functions Gn(x,y,t) and G A (x,y,t) are 
related to the transition amplitude by 

G R (x,y,t) = e(t)U(x,y,t), 

G A (x,y,t) = -9{-t)U{x,y,t), (134) 

and these satisfy the inhomogeneous equation 

(ihd t -Hx)G (x, y, t) = ihd (t) 5 3 (x) . (135) 

If we can solve this equation then we can immediately construct the transition am- 
plitude using the relation 

U (x, y, t) = G R (x, y, t) - G A (x, y, t) . (136) 

9.1 Example: the free Newtonian particle 

Given the Hamiltonian 

H = ™ < 137 > 

in continuous time mechanics we can readily find the Greens functions in the quantum 
theory. We find 

^ . , „ . . / — im \ 3 / 2 f i m (x — y) . (x — y) 1 

G « (x - y - () = »W(5ri«) exp {ft V - }' 

G.fx.y.i) = -«(-f)y expj- »' y ' |. (138) 

from which we construct the transition amplitude 

/— im\ 3 / 2 f i m (x — y) . (x — y) "1 

u ^ = (^m) Ha L ■ t>0 - (139) 



This satisfies the closure condition 

J d 3 yU (x, y; T) U* (z, y; T) = 5 3 (x - z) (140) 

for a system amplitude and demonstrates the essential point that it is possible to 
construct examples of discrete time quantum mechanics from continuous time quan- 
tum mechanics. The converse need not be true. Given a system amplitude, it may 
be impossible to find a compatible operator equivalent to some second order Hamil- 
tonian operator in continuous time mechanics. It is not difficult to find examples of 
normal systems where this occurs. 
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10 The discrete time harmonic oscillator 



Given the discrete time harmonic oscillator system function 

F n = F (x n+1 , x n ) = \a (x 2 n+1 + a£) - (3x n+1 x % 



(141) 



we note that is an example of a normal system. This leads us to define the system 
amplitude to be 



U n (x,y) = kexp(iF(x,y)/h) 



k exp 



2h 



ax 2 + ay 2 — 2(3xy 



(142) 



where k is some constant. From the unitarity condition ([97]) we find ( |1 1 2| ) and from 
( |108|) we find (|113|) , so we see that the above system co-ordinates are indeed normal. 
The self adjoint diagonal operator with operator component 



Cl = y- 1 \-h 2 d 2 x + (f3 2 - a 2 )x 2 



(143) 



is compatible with the system amplitude ( |142| ) . 

The interpretation of this is that C x is the operator corresponding to the Logan 
constant for the classical discrete time harmonic oscillator 



axy. 



(144) 



To see this explicitly, consider the operators x n and x n +i- The Logan constant for the 
harmonic oscillator is quantised according to the standard rule 



C (3 (x n x n -\- x 



(145) 



A suitable co-ordinate representation of these operators with respect to the basis B n 
is 

(146) 



Xn — >X. 



Xn+l — y V X ~ ^~£p x 



and then the operator (|145|) is represented by (|143| ) . 

We see from the potential term in the differential operator ( |143[ ) that a complete 
set of physical states can be found as eigenstates of the operator provided (3 2 > a 2 . 
This corresponds precisely to the elliptic region discussed in the classical theory. 

If the constants satisfy the elliptic condition then we may construct annihilation 
and creation operators for the system. These are diagonal with respect to any of the 
bases B n and are given by 



in 9 



Xn+1 ^ 2* 



And 



dx \x, n) 



h 



l-rfx+ -d x \ (x,n\ 



a' = 



-xe 



-iO 



n 



dx \x, n) | yl — rfx — —d x 



(x, n \. 



(147) 
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These operators satisfy the commutation relation 



2hJl - V 2 , x 

= (148) 

Using the evolution relation ( |106| ) and the operator equation of motion 

x n+1 = 2r]x n - (149) 

we find 



(150) 



but this does not mean this operator is conserved. A conserved operator according 
to our definition must be compatible with the timestep operator U n . We find that the 
creation and annihilation operators satisfy the relations 

U n a n - e ie a n U n = 0, 
U n a+ - e- ie atU n = 0, (151) 

which is reminiscent of various deformed commutators encountered in q-mechanics. 
However, the operator ( |145| ) corresponding to the Logan constant C n = ^(3a* n a n is 
compatible with the timestep operator and is therefore an invariant of the motion. 
We find 



C = j(3{a + a + aa + } = i/3d + a + ^Jl — rfTi 



i 



{_^2 q2 1 /" > 

2 p + |/3 (l - ?? 2 ) x 2 \ (x,n\ = J dx \x,n)C x (x,n\. 

(152) 

This Logan invariant is a close analogue of the oscillator Hamiltonian in continuous 
time mechanics and the eigenstates of the former follow the same pattern as the 
eigenstates of the latter. For example, there is a ground state |\l/o) satisfying the 
relation 

a n |*o) = (153) 

with normalisable wave function \I/ (x) = ^ (0) exp j— jfly/l — r] 2 x 2 /h} . This wave 
function is also an eigenstate of the Logan invariant operator, with 



C%(x) = yi-r 1 2 h^o(x). (154) 

These results hold only for |^| < 1. We note from §5.3 that f3 = m(6 + T 2 u 2 )/6T 
and that 

c h 2 

lim — = d 2 + -muj 2 (155) 

when we identify our system with the Newtonian oscillator. 
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11 The inhomogeneous oscillator 



In this section we discuss the inhomogeneous harmonic oscillator, which serves as 
a prototype for the application of our quantisation principles to field theories. We 
will use the source functional techniques of Schwinger to obtain the ground state 
functional and various n-point functions of interest. Because the Schwinger method 
deals with time ordered products we should expect the discretisation of the time 
parameter to involve some changes in the details of the calculations. 

First, given a system function F n = F (x n , x n+1 ) we are free to introduce an 
external source in any convenient way, as ultimately this will be set to zero. Our 
choice is to define the system function F n [j] in the presence of the external source as 

F n = F n + ^ T j n+1 x n+1 + \T Jn x n) (156) 

a choice which allows the construction of time ordered product expectation values 
directly. The action sum from time MT to time NT (N > M) now becomes 

A NM {j] = A NM + I TjmXm + l - T]N X N + T JnX n , M < N (157) 

n=M+l 

from which Cadzow's equation of motion is found to be 

d 



dx r . 



| F n + F n-i j + Jn = o, M < n < N. (158) 



Quantisation is introduced via the Schwinger action principle modified for discrete 
time. We postulate that for an infinitesimal variation 8A NM [j] of the action operator 
then 

5(0, N\if>, M) j = U<f),N\5A NM \j\ \*p,M) j , M < N (159) 

for any states \(f),N), \ip,M) at times NT, MT respectively, with evolution in the 
presence of the source. Independent variation of the j n for M < n < N then leads to 
the equations 

^ — (^N^MY = '-(^NlxM^My 

— til f) 

-— — (<j>,N\iP,M¥ = (<j>,N\x n \il>,My, M < n < N 

" " (<f>,N\4>,My = U<t>,N\x N \il),My. (160) 



T dj N 



Further application of the principle leads to expectation values of time ordered prod- 
uct of operators, such as 



-ih\ 2 d 

-jrj Qj^ih N\ip, MY = (<!>, N\Tx m x n \iP, MY, M<m,n<N (161) 



)2 
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where the symbol T denotes discrete time ordering. For example, 

Tx m x n I Q m — n -\- —5 m — n I x m x n -\- ( & n — m -\- —5 m — n I x n x 



O f'" ■ — n I Tn^ri, 1 \ — it — III 1 2 



where n is the discrete step function, defined by 



6 n = +1, n > 0, 

=0, n < (163) 



and 5 n is the Kronecker delta, defined by 



+1, n = 0, 

0, n ^ 0. (164) 



Given the harmonic oscillator system function 

f 

2T 6 



pn _ m (xn+i - x n f Tmu 2 , 2 2 \ 

F cyrp r- \ X n+l + %n+lXn + % n J 



(165) 



then the classical discrete time harmonic oscillator in the presence of the external 
source j n satisfies the equation 



where r] = a//3 with 



T 

x n+ i = 2r]x n - + — jn, (166) 

c P 



6T 6T V ; 



As discussed previously, elliptic (oscillatory) solutions occur for rf < 1 whereas hy- 
perbolic solutions occur for rf > 1. We will now discuss these possibilities separately. 

11.1 The elliptic regime 

The importance of the elliptic regime rj 2 < 1 stems from the fact that in field theory 
this corresponds to physical particle configurations of the fields, i.e., solutions which 
can be normalised. 

Now define the action of the {classical) discrete time displacement operator U n by 
the rule 

U n f n = f n+ i (168) 



for any function of the index n, where n is real. Then ( |166| ) may be written in the 
form 

(U n -2r ] + U- 1 )x n = ^ 3n . (169) 
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To solve (|166)) for the elliptic case we first define the following: since rf < 1 we 
write c = cos(#) = rj and s = sin (8) = — r] 2 > 0, taking < 6 < ir. If we 

define s a = sin (a9) where a is real then a useful identity is 



From this we deduce 
which is equivalent to 
We define the matrices 

A™ : 

and use ( |170| ) to prove 
If we write 



s a-\-l + s a-l = 2c5 a , 

(u a -2 V + u- 1 ) s a = o. 



A = A 1 , 



T ■ 

ajn 



then (|166|) may be written in the form 

X n = AX n _i + 3 n . 



(170) 

(171) 
(172) 

(173) 

(174) 
(175) 

(176) 



This equation may be readily solved using the properties of the s a functions and by 
diagonalising the matrix A. We choose Feynman boundary conditions, specifying the 
particle to be at position i M in the past (at time MT) and at position in the 
future (at time NT), giving 



SSN-M^n — SSN-n%M + ss n-M x N 
n-l 



T 

1 



N 



s N~n s M-mjm + s N-n s M-njn + ^ s M-n s N-mjm ( ,(177) 
m=M m=l+n 



which is valid only for M < n < N . This can be tidied up into the form 

N 

_ sn-uXm sm-uXn _ T sr^ r< nm o 



where 



r~mm 



SN-M sm-n 



s N-n s M-m 



;i78) 



Pss f 



N—M 
SN-nSM-m 

@ SS n-m 
s M-n s N—m 



(3ssn- 



m=M 

M < m < n < N 
M < m = n < N 
M < n < m < N. 



(179) 



M 
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Then Gjvm satisfies the inhomogeneous equation 

P (U n -2 V + U- 1 ) G n N m M = -5 n . m , M < n < N. (180) 

Up to this stage we have taken N > M with both finite, but normally we will 
be interested in the scattering limit N — > +00, M — ► —00. Also, we have appeared 
to have overlooked the possibility that s N _ M vanishes in the denominator of the 
propagator ( |179| ) for some values of N and M. We shall address both of these issues 
directly now. 

Our method of avoiding possible singularities is to extend the Feynman — ie pre- 
scription to the 9 parameter. By inspection of the equation 

^™0) = -6TTVT ( 181 ) 

we deduce that 

u 2 -> uj 2 - it 6 -> - it, rj -> 77 + ie. (182) 

With this deformation of the # parameter and taking the limit iV — > +00, M — » —00, 
we find 

00 

x n = x n -T X! C™- m J m (183) 

m=— 00 

where x„ satisfies the homogeneous equation 

(C7 n - 2 V + U" 1 ) x n = (184) 

and 

G n f m = (e i(m -" )e 9 n _ m + 5 m _ n + e j( "- m)e 6 m _ n ) . (185) 

This is the discrete time analogue of the harmonic oscillator Feynman propagator 
and reduces to it in the limit T — > 0, nT — ► t. A direct application of the discrete 
time Schwinger action principle to the operator equation of motion 

(U n -2 V + U- 1 ) x n = ^ Jn (186) 



then gives the ground state vacuum functional 



Z [7] = Z [0] exp — - £ J n G n F - m 3m \, (187) 

[ n,m=-oo J 

essentially solving the quantum problem. Using this result and ( |16U| ) we find that in 
the limit j —* 

(0\x n x n \0) = A (0|x n+1 x n |0) = ^e"". (188) 
From this we deduce 

(0\[x n+1 ,x n )\0) = ^-, (189) 
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which agrees exactly with the discrete time oscillator commutation relation 



[x n+1 ,x n ] = —j^- (190) 

found previously. 

Further ground state expectation values of commutators may be obtained by using 
the result 

(0|Tx m x„|0) = %hG°f m = ^— a e' l[n - Me . (191) 

Zp sin v 

For example, we find 

(0|[x n+2 ,x n ]|0) = ^^, (192) 
which agrees with the commutator 

[x n+2 ,x n } = — - — (193) 

obtained from the operator equation of motion (|149|) and the commutator (|190|) . 
Another verification of the consistency of our methods is that we may use (|188 ) 



directly to find the ground state expectation value of the Logan invariant ( |145| ) for 
the discrete time harmonic oscillator. We find 



(0|(7"|0) = ^ v /l-r / 2, (194) 
which agrees exactly with previous results. 

11.2 The hyperbolic regime 



Because Tu is real and positive the controlling parameter rj as given by ( |181|) takes 
vales only in the regions 



elliptic 
parabolic 
hyperbolic 



- 1 < j] < 1 : 0<Tu < 2V"3 

77 = —1 : Too = 2V3 
-oo<77<-1: 2 v / 3<Tw. (195) 



If we parametrise 77 by the rule 77 = cos (z) where z is complex then if we take 

T) = cos 6: < Tuj < 2y/3 

7] = - cosh A : 2\/3 < Tuj, (196) 



then the range of possibilities (|195|) corresponds to a contour T in the complex z 



8 — iX plane which runs just below the real axis from the origin to ir and then runs 
from 71 to 71 — zoo. The elliptic regime corresponds to values of z on the first part of 
the contour, for which A = + e, where e is infinitesimal and positive, corresponding 
to the Feynman — ie prescription. 
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The hyperbolic region corresponds to the part of the contour given by z = n — 
iX, : A > 0. For this region analytic continuation of the s n functions leads to 

s n ^i{-l) n+1 ~ Sn (197) 

where s n = sinh (nA) . From this the analytic continuation of the finite interval prop- 
agator (|179|) gives 

/ \n—m 

G 'nA'I = n — {s N —nS M -rnS) n—m ~\~ §N—m§M—m,3n—m ~\~ Qm-nS M -n$ N -m} > 

M < n,m < N, (198) 

which satisfies the equation 

(3 {U n - 2?? + U- 1 }GTm = Sn-m- (199) 
Taking the limit iV = — M — > oo gives the infinite interval propagator 

/ \l+n— m 

G n f m = l ~ J {e^ x Q n . m + 5 n . m + e^- m ) A 9 m _ n } (200) 

which satisfies the equation 

P [U n -2 V + U- 1 } G n f m = -5 n . m . (201) 

11.3 Comment 

The elliptic and hyperbolic Feynman boundary condition propagators can be sum- 
marised in the analytic form 

T(2 + cos z) ^_ iHz 



bmt sin z 
T (2 + cos z) 



6mi sin z 



(e~ inz e n + 5 n + e m2 6_ n ) , (202) 



where r] = cos z and z lies somewhere on the contour V discussed above. Then A^ (77) 
satisfies the equation 

P {U n -2 V + U- 1 } A™ ( V ) = -5 n . (203) 

However, physical states correspond only to points in the elliptic regime, as the 
wave functions will not remain normalisable in time otherwise. For example, we 
showed previously that the ground state wavefunction for the discrete time oscillator 
in the elliptic regime is given by 



^0 0) = ^0 (0) exp ~pJl - rfx 2 /h (204) 



demonstrating that analytic continuation to the hyperbolic regime will not give a 
normalisable ground state wave-function. 
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12 Concluding remarks 



In this paper we have shown that a consistent approach to the discretisation of time 
results in a consistent dynamical framework. Once the initial psychological hurdle 
of accepting a dynamics without time derivatives has been jumped, then such a 
theory becomes as reasonable as continuous time dynamics. Indeed, by taking the 
fundamental time interval T small enough, it would appear possible to duplicate 
or approximate conventional theory as closely as required. The obvious question, 
why consider discrete time mechanics at all? has two replies. First, it may be 
the case after all that there is some fundamental limit to time intervals, and so it 
becomes a matter of curiosity as to how far we can go along that road. Second, there 
may be some novel properties in this approach which could prevent or alleviate the 
notorious problems with divergences which plague conventional field theories. It is 
worth considering any approach to the regularisation in field theory which is based 
on just one assumption, namely that of a discrete time. The behaviour of the discrete 
time oscillator holds the promise of potentially useful properties which may provide 
a cutoff for particle energy. 

In the next paper we apply our methods to classical field theories. However 
no divergence problems appear at that stage. In the third paper of this series, we 
shall consider second quantisation, and the issues to do with divergences of Feynman 
diagrams will be discussed in subsequent papers in some detail. 
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